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Abstract 



We propose a new analytical method to solve for the nonexactly 

"^ ' solvable Schrodinger equation. Successfully, it is applied to a class of 

Qj. spiked harmonic oscillators and truncated Coulomb potentials. The 

"i I utility of this method could be extended to study other systems of 

-y ■ atomic, molecular and nuclear physics interest. 
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In atomic, molecular and nuclear physics, spiked harmonic oscillators and 
truncated Coulomb potentials are of significant interest. Realistic interaction 
potentials often have a usually repulsive core [1-5]. The simplest model of 
such a core is provided by the spiked harmonic oscillators 

V{q) = Ciq^ + C2q-\ Ci,C2,6>0. (1) 

On the other hand, the truncated Coulomb potential has been founded to 
be pertinent in the study of the energy levels of the hydrogen - like atoms 
exposed to intense laser radiation [6-11]. It has been shown [8,10] that under 
Kramers - Henneberger transformation [12] the laser - dressed binding poten- 
tial for the hydrogenic system, often called laser - dressed Coulomb potential, 
may be well simulated by 

Where the truncation parameter c is related to the strength of the irradiating 
laser field. 

Thus it is interesting to carry out systematic studies of the bound - states 
of these potentials. Hall and Saad [4] have studied the spiked harmonic os- 
cillator potentials via smooth transformations method (STM) of the exactly 
solvable potential V{q) = C2q'^ + C2g~^ to obtain lower and/or upper energy 
bounds. They have also calculated the energy eigenvalues using direct nu- 
merical integrations of Schrodinger equation [4]. Dutt et al [6] have used a 
shifted 1/N expansion technique (SLNT) to carry out the energy levels of the 
laser - dressed Coulomb potential and compared their results also with those 
of direct numerical integrations [7]. Nevertheless, neither SLNT nor STM is 
utilitarian in terms of calculating the eigenvalues and eigenfunctions in one 
batch. Because of the complexity in handling large - order corrections of the 
standard Rayleigh - Schrodinger perturbation theory, only low - order cal- 
culations have been reported for SLNT [6,13] and large - order calculations 
have been neglected. Eventually, the results of SLNT are not as accurate as 
sought after. 

In this paper we formulate a method for solving Schrodinger equation. In 
one batch, one should be able to study not only the eigenvalues but also the 
eigenfunctions. It simply consists of using 1// as a perturbation expansion 
parameter. Where, / = / — /3, / is a quantum number, and /3 is a suitable 



shift introduced to avoid the trivial case / = 0. Hence, hereinafter, it should 
be called pseudoperturbative shifted - / expansion technique (PSLET). 

The construction of our method starts with the time - independent one - 
dimensional form of Schrodinger equation, in h = m = 1 units. 
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Where the quantum number / may specify parity, (—1)'+^, in one - dimension 
( l=-l or /=0, and q G (— oo, oo)) or angular momentum in three - dimensions 
( /=0, 1, ■ ■ ■ , and q G (0, oo)), and nr=0, 1, ■ ■ ■ counts the nodal zeros [6,14- 
17]. 

To avoid the trivial case /=0, the quantum number / is shifted through 
the relation I = I — (3. Eq.(3) thus becomes 

\l'^ + Viq)\ ^„.,Kg) = i?„.,i^„.,/(g), (4) 

vM^ ^"'^^^g"^'^"^' .^^M. (5) 

Herein, it should be noted that Q is a constant that scales the potential V{q) 
at large - / limit and is set, for any specific choice of / and Ur, equal to /^ at 
the end of the calculations [13,14]. And, f3 is to be determined in the sequel. 
Our systematic procedure begins with shifting the origin of the coordi- 
nate through 

x = l'/'{q-qo)/qo, (6) 

where qo is currently an arbitrary point to perform Taylor expansions about, 
with its particular value to be determined. Expansions about this point, 
a; = (i.e. q = qo), yield 

y n=0 Vo 
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Obviously, the expansions in (7) and (8) center the problem at an arbitrary 
point qo and the derivatives, in effect, contain information not only at qo but 



also at any point on the axis, in accordance with Taylor's theorem. Also it 
should be mentioned here that the scaled coordinate, equation (6), has no 
effect on the energy eigenvalues, which are coordinate - independent. It just 
facilitates the calculations of both the energy eigenvalues and eigenf unctions. 
It is also convenient to expand E as 



En.,l = E Kr- 



n=-2 



Equation (4) thus becomes 
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where the prime of V{qo) denotes derivative with respect to qo- Equation 
(10) is exactly of the type of Schrodinger equation for one - dimensional an- 
harmonic oscillator 
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where P{x) is a perturbation - like term and Eq is a constant. A simple 
comparison between Eqs.(lO), (11) and (12) implies 
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Equations (14) and (15) yield 
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Here Qo is chosen to minimize E!^^ /, i. e. 



dE. 



(-2) 
rir.l 



'2 f(-2) 



anrf 



c^^£;. 



Urd 



>0. 



(15) 

(16) 

(17) 

(18) 
(19) 

(20) 



Hereby, V{q) is assumed to be well behaved so that E^ 2) j-^as a minimum go 
and there are well - defined bound - states. Equation (20) in turn gives, with 



l-P=^qlV'{qo). (21) 

Consequently, the second term in Eq.(ll) vanishes and the first term adds 
a constant to the energy eigenvalues. It should be noted that energy term 
PE^ I has its counterpart in classical mechanics. It corresponds roughly to 
the energy of a classical particle with angular momentum Lz=l executing 
circular motion of radius Qo in the potential V{qo). This term thus identifies 



the leading - order approximation, to all eigenvalues, as a classical approx- 
imation and the higher - order corrections as quantum fluctuations around 
the minimum g^, organized in inverse powers of /. 

The next leading correction to the energy series, lE~i, consists of a 
constant term and the exact eigenvalues of the unperturbed harmonic oscil- 
lator potential w'^x'^/2. The shifting parameter (3 is determined by choos- 
ing JeI^ / =0. This choice is physically motivated. It requires not only the 
agreements between PSLET eigenvalues and the exact known ones for the 
harmonic oscillator and Coulomb potentials but also between the eigenfunc- 
tions as well. Hence 



P = - 



- + K + -)w 



where 



w 



3 + 



qoV'iqo) 



\ V'{q^) ■ 



Then equation (11) reduces to 
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Equation (10) thus becomes 
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When setting the nodeless, n^ = 0, wave functions as 

^^oMq)) = exp{Uo,i{x)), (29) 

equation (28) is readily transformed into the following Riccati equation: 

- l[u"{x) + u\x)u'ix)] + £ .W(x)r"/^ = i (^(^ + Ar) 
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Hereinafter, we shall use U{x) instead of Uo,i{x) for simplicity, and the prime 
of U{x) denotes derivative with respect to x. It is evident that this equation 
admits solution of the form 
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Substituting equations (31) - (33) into equation (30) implies 
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where primes of U^'^\x) and G*^"'''(a;) denote derivatives with respect to x. 
Equating the coefficients of the same powers of / and x, respectively, ( of 
course the other way around would work equally well) one obtains 
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and integration over dx yields 
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■ ■ ■ and so on. Thus, one can calculate the energy eigenvalue and the eigen- 
functions from the knowledge of Cm,n and Dm,n in a hierarchical manner. 
Nevertheless, the procedure just described is suitable for systematic calcula- 
tions using software packages (such as MATHEMATICA, MAPLE, or RE- 
DUCE) to determine the energy eigenvalue and eigenfunction corrections up 
to any order of the pseudoperturbation series. 

Although the energy series, Eq.(9), could appear divergent, or, at best, 
asymptotic for small /, one can still calculate the eigenenergies to a very good 
accuracy by forming the sophisticated Pade' approximation to the energy 
series. The energy series, Eq.(9), is calculated up to -Eq//P by 

Eo,i = PEJ^f + </ + ■ ■ ■ + <)/r^ + 0{1/F), (51) 

and with the P^{1/1) and P^{1/1) Pade' approximants it becomes 

i?o,^[3,3]=PEj;') + P|(l/r). (52) 

and 

EoASA] = PEir,'^ + p,\i/i). (53) 



Hereby, an " if " statement is in point. If the energy series, eq.(9), is a Stieltjes 
series, though it is difficult to prove, then -Eo,/[3, 3] and -Eo,«[3, 4] provide upper 
and lower bounds to the energy [18,19]. Our strategy is therefore clear. 
Let us begin with the spiked harmonic oscillators 

V{q) = \{q' + aq-') (54) 

for which Eq.(22), with rir = 0, implies 



1 8qo + ab{b - 2)qo 

In turn Eq.(21) reads 
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Equation (56) is explicit in qo and evidently a closed form solution for qo 
is hard to find, though almost impossible. However, numerical solutions are 
feasible. Once qo is determined the coefficients Cm,n and Dm^n are obtained in 
a sequel manner. Consequently, the eigenvalues, Eq.(51), and eigenfunctions, 
Eqs.(31)-(33), are calculated in the same batch for each value of a, b, and /. 
In tables 1 and 2 we list PSLET results Ep, Eq.(51), along with [3,3] and 
[3,4] Pade approximants, Eqs.(52) and (53) respectively. The results of the 
smooth transformations method (STM) [4] and direct numerical integration 
(DNI) [4] are also displayed for comparison purposes. 

Our calculated values of the bound - state energies, Ep, compare well 
with those from direct numerical integrations [4]. In table 1 the Fade ap- 
proximants -E[3,3] and -E[3,4] are almost in total agreement with those of 
Hall and Saad [4] from DNI of the Schrodinger equation. Moreover, it is 
evident that -E[3,3] and -E[3,4] have provided upper and lower bounds, re- 
spectively, to the energy series. However, the same can not be concluded 
from table 2. Eventually, our computed values of the bound - state energies, 
Ep, do not contradict with the upper and/or lower bounds reported by Hall 
and Saad [4] from the smooth transformations method (STM). 

Moreover, our result for b = 2 listed in table 1 is in excellent agreement 
with the exact one 65.2534584 obtained from Eq.(2) of ref. [4]. On the other 
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hand, one would rewrite the centrifugal term in (3) plus the potential (54) 
as /'(/' + l)/(2g2) + g2/2, where /' = -1/2 + ^{l + l/2f + a, and proceed by 
shifting the irrational quantum number /' through I = I' — j3. In this case, one 
obtains the known exact result Ep = (/' + 3/2) for the harmonic oscillator 
g^/2 from the leading term t^E^^"^' and the remainder energy corrections are 
identically zero. 

Next, we consider the laser - dressed Coulomb potential 



V{q) 



^/q^Tc^ 



c> 0. 



(57) 



In this case 
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Again, we numerically solve for g^, and proceed exactly as above to calculate 
the energy eigenvalues and eigenfunctions in the same batch. In tables 3 and 
4 we collect the results for the truncation parameter c = 1,5, 10, 50, 100, 200 
based on our approach. The energies Ep, Eq.(51), compare well with those 
of Singh et al. [7] from numerical integrations. The Pade approximants 
i?[3, 3] and i?[3, 4] are in almost complete accord with those of Singh et al.[7]. 
However, they do not provide upper and lower bounds to the energy series, 
Eq.(51). Perhaps, it should be mentioned that the approximate binding po- 
tential Eq.(57) is valid for a hydrogen atom in a laser field which corresponds 
to a truncation parameter c in the range 20-60 [6]. Higher and lower values 
of c have been considered for academic interest only. 

Before we conclude some remarks deserve to be mentioned. 

For the two problems discussed in this paper, we have shown that it is an 
easy task to implement PSLET without having to worry about the ranges of 
couplings and forms of perturbations in the potential involved. In contrast to 
the textbook Rayleigh - Schrodinger perturbation theory, an easy feasibility 
of computation of the eigenvalues and eigenfunctions, in one batch, has been 
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demonstrated, and satisfactory accuracies have been obtained. Moreover, a 
nice numerical trend of convergence has been achieved. Nevertheless, another 
suitable criterion for choosing the value of the shift /9, reported in Ref. [14], 
is also feasible. This reference should be consulted for more details. 

It is not easy to prove that the energy series Eq.(51) is a Stieltjes series. 
But, if it is a Stieltjes series, the [N, N] and [N, N + 1] Fade approximants 
provide upper and lower bounds to the energy series. Table 1 bears this 
out. Moreover, in view of the results listed in tables 1-4 one can confidently 
conclude that the [3,3] and [3,4] Fade approximants to the energy series 
Eq.(51) can be used to determine the energy eigenvalues to a very satisfactory 
accuracy. 

From the knowledge of Cm,n and Dm,n one can calculate, in the same 
batch, the wave functions to study electronic transitions and multiphoton 
emission occurring in atomic systems in the presence of intense laser fields, for 
example. Such studies already lie beyond the scope of our present methodical 
proposal. 

Finally, the attendant technique FSLET could be applied to Schrodinger 
equation with rational potentials, such as the nonpolynomial oscillator V{q) = 
q^ + Ag^/(1 + gq^). This type of potential is an interesting model in laser and 
quantum field theories [20]. The feasibility of FSLET extends also to a class 
of screened Coulomb potentials, which have relevance in atomic and plasma 
physics, and to some other models of interest [21-26, and references therein]. 
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Table 1: Is - state energies, in % 



m 



1 units, of the potential V{q) 



{q^ + 1000/g'')/2. Where Ep represents PSLET results, Eq.(51), E^'^ with 
U and L denote upper and lower bounds from STM [4] , and En from DNI 
[4]. -E[3,4] is the [3,4] Pade approximant obtained by replacing the last j 
digits oi E[3, 3] with the j digits in parentheses. 



b 


Ep 


E3,3] & (^3,4]) 


Es 


En 


0.5 


415.88978 


415.889786 (86) 


416.30977^ 


415.88979 


1.0 


190.72330 


190.723308 (07) 


190.99213^ 


190.72331 


1.5 


104.41022 


104.410224 (24) 


104.53993^ 


104.41022 


1.9 


71.06157 


71.0615789 (87) 


71.08686^ 


71.06158 


2.0 


65.25345 


65.2534589 (86) 


65.25346 


65.25346 


2.1 


60.15200 


60.1520114 (11) 


60.12704^ 


60.15201 


2.5 


44.95547 


44.9554855 (50) 


44.83349^ 


44.95549 


3.0 


33.31675 


33.3167621 (18) 


33.07940^ 


33.31676 


3.5 


26.10884 


26.1088462 (48) 


25.76204^ 


26.10885 


4.0 


21.36950 


21.3694640 (14) 


20.91865^ 


21.36964 


4.5 


18.10194 


18.1018377 (10) 


17.55218^ 


18.10183 


5.0 


15.76134 


15.761144 (25) 


15.11758^ 


15.76113 


5.5 


14.03138 


14.03112 (07) 


13.29842^ 


14.03107 


6.0 


12.71886 


12.71879 (61) 


11.90153^ 


12.71862 
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Table 2: Is - state energies, in h = m = 1 units, of the potential V{q) = 
{q^ + a/g^/^)/2. Where Ep represents PSLET results, Eq.(51), Es denotes 
the lower bounds from STM [4] , and E^ from DNI [4]. ^[3,4] is the [3,4] 
Pade approximant obtained by replacing the last j digits of E[3, 3] with the 
j digits in parentheses. 



a 


Ep 


E3,3 L{E3,4) 


Es 


En 


1000 


44.95547 


44.9554855 (50) 


44.83349 


44.95549 


100 


17.54168 


17.541911 (899) 


17.41900 


17.54189 


10 


7.73423 


7.73606 (548) 


7.61169 


7.73511 


5 


6.29679 


6.29988 (756) 


6.17394 


6.29647 


1 


4.32861 


4.528 (290) 


4.20453 


4.31731 


0.5 


3.85740 


3.8308 (289) 


3.74611 


3.84855 


0.05 


3.13431 


3.1606 (893) 


3.10954 


3.15243 


0.005 


3.01445 


3.0199 (201) 


3.01178 


3.01905 
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Table 3: Bound - state energies, in h = m = 1 units, of the potential 
V{q) = — (g^ + c^)^^/^ for the Is, 2p, 3d, and 4f states. Where Ep represents 
PSLET results Eq.(51), Eslnt from SLNT [6], and En from DNI [7]. E[3,A] 
is the [3,4] Pade approximant obtained by replacing the last j digits oi E[3, 3] 
with the j digits in parentheses. 



c State 


-Ep 


-^[3,3] & {-E[3,4]) 


—Eslnt 


—En 


1 Is 


0.27412 


0.27478 (62) 


0.27596 


0.27439 


2p 


0.113087 


0.11296 (303) 


0.112826 


0.113024 


3d 


0.0544357 


0.0544371 (82) 


0.054442 


0.0544362 


4f 


0.03106845 


0.03106846 (47) 


0.031069 


0.03106846 


5 Is 


0.1070836 


0.1070813 (10) 


0.107396 


0.1070814 


2p 


0.06819140 


0.06818667 (33) 


0.068233 


0.06818716 


3d 


0.04325586 


0.04325730 (20) 


0.043247 


0.04325755 


4f 


0.02810534 


0.02810520 (25) 


0.028101 


0.02810524 


10 Is 


0.06373831 


0.06373817 (21) 


0.063820 


0.0637389 


2p 


0.04620043 


0.04619903 (00) 


0.046228 


0.04619904 


3d 


0.03315868 


0.03315855 (53) 


0.033164 


0.03315859 


4f 


0.02380662 


0.02380672 (71) 


0.023806 


0.02380674 
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Table 4: Bound - state energies, in h = m = 1 units, of the potential 
V{q) = — (g^ + c^)^^/^ for the Is, 2p, 3d, and 4f states. Where Ep represents 
PSLET results Eq.(51), Eslnt from SLNT [6], and En from DNI [7]. E[3,4] 
is the [3,4] Pade approximant obtained by replacing the last j digits oi E[3, 3] 
with the j digits in parentheses. 



c 


State 


~Ep 


-E[3,3] k (-_ 


E[3,4]) 


—Eslnt 


—En 


50 


Is 


0.01626071 


0.0.01626072 


(71) 


0.016263 


0.01626072 




2p 


0.01408837 


0.01408837 


(37) 


0.014090 


0.01408838 




3d 


0.01215871 


0.01215871 


(71) 


0.012160 


0.01215871 




4f 


0.01045842 


0.01045842 


(42) 


0.010459 


0.01045842 


100 


Is 


0.00862978 


0.00862978 


(78) 


0.008630 


0.00862978 




2p 


0.00780013 


0.00780013 


(13) 


0.007800 


0.00780013 




3d 


0.00703519 


0.00703519 


(19) 


0.007035 


0.00703519 




4f 


0.00633273 


0.00633273 


(73) 


0.006333 


0.00633273 


200 


Is 


0.00450285 


0.00450285 


(85) 


0.004503 


0.00450286 




2p 


0.00419307 


0.00419307 


(07) 


0.004193 


0.00419307 




3d 


0.00390020 


0.00390020 


(20) 


0.003900 


0.00390020 




4f 


0.00362385 


0.00362385 


(85) 


0.003624 


0.00362385 



